Two-sphere low Reynold's propeller 
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A three-dimensional model of a low-Reynold's swimmer is introduced and analyzed in this paper. 
This model consists of two large and small spheres connected by two perpendicular thin rods. The 
geometry of this system is motivated by the microorganisms that use a single tail to swim, the 
large sphere represents the head of microorganism and the small sphere resembles its tail. Each rod 
changes its length and orientation in a non-reciprocal manner that effectively propel the system. 
Translational and rotational velocities of the swimmer are studied for different values of parameters. 
Our findings show that by changing the parameters we can adjust both the velocity and the direction 
of motion of the swimmer. 

PACS numbers: 87.19.ru, 47.15. G-, 45.40.Ln 



I. INTRODUCTION 

The propulsive motion of artificial and biological 
micron-scale objects is an interesting problem at low 
Reynolds hydrodynamics. In this condition the dynam- 
ics is dominated by viscous forces. Examples of these 
micron scale objects include biological microorganisms 
like bacteria and also man-made micro-swimmers, useful 
to operate at microfiuidic investigations [l[. 

Propulsive motion at low Reynold's number is subject 
to the Scallop theorem 0]. At small scales, where, the 
Reynolds number is very low, the governing hydrody- 
namic equations, i.e. the Stokes and continuity equa- 
tions are linear and invariant under time reversal Q. 
Any reciprocal shape deformation retraces its trajectory 
and the system stay back at the point where it started. 
In order to achieve a net translational displacement, 
the system should perform the body deformations in a 
non-reciprocal manner. As mentioned by Purcell a low- 
Reynolds propeller must have at least two internal de- 
grees of freedom and he proposed a three-link swimmer. 
The detailed motion of Purcell's swimmer was examined 
by H. Stone where, it was shown that Purcell's system 
could swim and its dynamical properties were calculated 
Inspired by Purcell's system, a low Reynolds swim- 
mer constructed by three linked spheres was introduced 
and analyzed by Najafi et al. @ , and experimentally real- 
ized by M. Leoni et al. @ . After Purcell's proposal there 
have been considerable scientific efforts in designing arti- 
ficial swimmers. Such swimmers would be useful in devel- 
oping microfiuidic experiments. Furthermore, progresses 
in assembling micro-swimmers show the possibility of us- 
ing micro-machines inside the biological cells for non in- 
vasive therapeutic treatments Q. On the other hand, 
there are many theoretical works devoted to the study 
of different aspects in the motion of biological microor- 
ganisms at low Reynolds number condition (8rJl2l]. Such 



FIG. 1: Schematic showing the geometry of a two-sphere 
swimmer. Two large and small spheres are connected through 
two perpendicular rods, one with fixed but the other with vari- 
able length. The short rod is rotating around the long rod. 
This model system resembles the motion of a bacterium that 
has a single tail. 



interests include sperm swimming, metachronal waves in 
cilia, E-Coli chemotaxis and coupling mediated by hy- 
drodynamic interaction between nearby microorganisms 
[HEl. For a review of recent progress on low Reynolds 
hydrodynamics of microor gan isms, see for example the 
review article by E. Lauga 15 1. 

Our first aim in this article is to present a simplified 
model that captures the characteristics of a swimming bi- 
ological organism like a bacterium. Dipolar far velocity 
field and asymmetric shape, corresponding to the head- 
tail geometry of the organisms are two important features 
of micro swimmers. We model these systems by consid- 
ering two spheres with different radii that are changing 
their separation. We will study the translational and an- 
gular motion of this system. 
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II. TWO-SPHERE MODEL 

Figure [T] shows the schematic geometry of a model 
swimmer composed of two spheres. As shown in this 
figure two small and large spheres with radii a and R are 
connected by two perpendicular and negligible diameter 
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FIG. 2: Trajectory of the small sphere seen in the frame of 
reference that is co-moving with large sphere. Here we chose 
h = 1 and ^ = 2. 



rods. Let denote the lengths of long and short rods by 
L and I respectively. The connection is established in a 
way that the angle between two rods is fixed to ^ while 
the relative angular position of small rod with respect 
to the large sphere can be varied. Additionally, we as- 
sume that the length of the long rod can be dynamically 
changed. In this case, the system will have two internal 
degrees of freedom, the length of the long rod L(t) and 
the rotational angle of the short rod <fi(t). 

The geometry which we are introducing here resem- 
bles the body shape of a bacterium with a single flag- 
ellum or cilium. Bacteria use beating patterns in their 
tails to move. The small sphere in our two-sphere model 
acts as a beating tail and the large sphere resembles the 
head of animal. The minimum condition for swimming 
at low Reynolds number can be achieved in our three- 
dimensional model. By changing the length of long rod 
and the angle of small one in a prescribed form we are 
able to choose the motion which breaks the time-reversal 
symmetry, the necessary condition for translational mo- 
tion and consequently propel the system. 

As an example for the internal motion of the system, 
we let the angle <fi(t) increase with constant angular veloc- 
ity and the length of long rod change periodically around 
an average length. The explicit form of this motion is 
given by: L(t) = L Q + /j, cos(w L i — O ), and <j)(t) = u^t. 
In this case, the position vector of the small sphere, seen 
in the reference frame that is co-moving with the large 
sphere, is given by: 

X = (I cos <j>(t), I sin <t>(t),L + h(t)), (1) 

where I is the length of the small rod and Lq represent 
the average length of the long rod. Figured] shows a typ- 
ical real space trajectory of the small sphere that is seen 
in the reference frame co-moving with the large sphere. 
Different choices for ujl and correspond to different 

forms of the internal motion. For — ^ = m = - with » and 

ui L q 1 

q as two positive integer numbers, we see that the paths 
in the (h, </>)-space are closed loops. One should note that 
the phase space of internal motion {{h, ^>)-space) is the 
surface of a cylinder. The axial direction on the cylin- 
der represents the /i-direction and the azimuthal angle is 
shown by the transverse direction on the cylinder. For 



m < 1, the phase space trajectory is a closed curve which 
traces exactly one complete turn round the cylinder while 
for m > 1 the trajectory is a closed loop that turns many 
times round the cylinder. In both cases, the geometry 
of the closed curves in the cylindrical-shape phase space 
are examples of nonreciprocal motions that can generate 
a net translational motion. 



III. POINT FORCE NEAR A RIGID SPHERE 

At zero Reynold's number the Stokes equations gov- 
ern the dynamics of fluid. The solution of the Stokes 
equation for a point force singularity is formulated in 
terms of Green's function and is called Stokcslct. For a 
point force singularity with strength f located at point 
X , the velocity field generated in the fluid is given by: 
u(X) = — ^G(X, Xo) • f , where rj is the viscosity of the 
fluid. C.W. Oseen has derived the explicit form of the 
Green's function G for an infinite flow that is bounded 
internally by a solid sphere with radius R For the 
explicit form of this solution, we refer to the original ar- 
ticle by Oseen. As it is manifested by Oseen's solution, 
the flow field due to a point force in the presence of a 
solid sphere can be regarded as the flow of the original 
point force and the flow due to the singular parts that 
are located at an image position inside the sphere. The 
location of the image point inside the sphere is given by 
its position vector, Xq = -^-Xo, relative to the sphere's 
center. 

As argued by Higdon the total force acting by a point 
force on fluid bounded by a no-slip sphere is equivalent to 
the total Stokeslet strength [I?} ■ Total Stokeslet strength 
includes the image point force inside the sphere. For a 
point force f , the image point force is defined by: f/ = 
(cyf r + c t f t ) where the radial and tangential components 
of this image force are given by: f r = f • Xo/Ao, ft = 
f — f r . Here two coefficients c r and Ct are given by: 

3 R 1 R 3 3 R 1 R 3 

In the next section we will use these results for analyzing 
the motion of two-sphere system. 

IV. TWO-SPHERE DYNAMICS 

In this section we will develop the dynamical equations 
for two-sphere system. To simplify the equations We will 
assume that the radius of small sphere is much smaller 
than the radius of large sphere (a <C R). We further 
assume that a is smaller than the characteristic distance 
between the spheres. With this approximation the veloc- 
ity field of the system is described by the velocity field 
of a point force that is moving near a rigid sphere. This 
simplification allows us to use the results of the preceding 
section and derive a more simpler dynamical equations 
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of the system. However one should note that the finite 
size effect of the small sphere can be systematically con- 
sidered by Faxen's theorem [l8j ]. 

To obtain the swimming velocity of the system we work 
in the reference frame, that is co-moving and rotating 
with the large sphere. In this coordinate system the ve- 
locity of the fluid at infinity is the swimming velocity. 
Denoting the swimming velocity of the system by V, and 
its angular velocity by 57, we can express the velocity field 
of the fluid at a general point X as: 

u(X) = -V-fi x X + M.V + fJ x m + G.f, (3) 

where the tensor M and vector m give the flow field due 
to the translational and rotational motion of a moving 
sphere. These quantities are given by: 

3 R rT XX . 1 i? 3 rT XX . i? 3 . 

In the absence of external force and torque, the swim- 
mer is force and torque free. Therefore we require that 
the total force and torque acting on the fluid by the sys- 
tem to be zero. Including the point force and its im- 
age counter part and adding the contributions due to the 
translational and rotational motion of large sphere we ar- 
rive at the following force and torque balance equations: 

J f + f/ + GmjRV = 

| X X f t - f^X X f t + 87T7/i? 3 ft = ^ 

The fluid velocity field at the location of small sphere 
is subject to the boundary condition: u(Xo) = Xq. To- 
gether with this boundary condition, the above equations 
make a complete set of dynamical governing equations for 
the swimmer. 

To solve the dynamical equations for the system, we 
can use the force and torque balance conditions and ob- 
tain a set of equations which relate the different com- 
ponents of the translational or angular velocities of the 
system to the component of the vector Xo in the follow- 
ing matrix form: 

V = A C 1 X 

n = b cr 1 x (6) 

where the detail for of the matrix elements [A]ij = a%j, 
[B]ij = bij and [C]ij = Cij are given at the appendix. 

V. RESULTS AND DISCUSSION 

In this section we will present the numerical solution 
for the governing equations and obtain the real space 
trajectory of the swimmer. For this purpose we plot 
the trajectory of large sphere. For the prescribed inter- 
nal motion given by (JlJ and a special choice of parame- 
ters, we have plotted figure El the space trajectory of the 
large sphere. As one can distinguish, the trajectory is a 
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FIG. 3: The trajectory of the swimmer in the (x, y, z)-space 
is plotted for two different values of (j>o- Other parameters are 
R = 1, a = 0.5, L = 4, ho = 0.1, I = 0.3 and w = 2ui = 1. 
Line show the real path of the large sphere. The swimmer 
starts its motion from the initial state where the large sphere 
is located at the origin and the long rod is orientated along 
the — z direction. As one can see, the overall swimming direc- 
tion can be varied by changing the parameters of the system. 
Average orientation of the long rod which is not shown here 
is along the average swimming direction. 
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FIG. 4: Average swimming velocity is plotted in terms of 
the length of long rod. Other parameters are set to: R = 1, 
a — 0.5, 4>o = 0, ho = 0.1, I = 0.3 and iv^ = ui = 1. 



helical-shaped path with overall translational movement 
in each round. The different characteristics of the tra- 
jectory, preferred direction, average swimming velocity 
and the effective radius of the helice can be controlled by 
the geometrical as well as dynamical parameters of the 
swimmer. 

The average orientation of the long rod, which is not 
shown in figure, is in the direction of the longitudinal 
axis of helix. This is achieved by numerically solving 
for the rotational velocity. Controlling and adjusting the 
dynamical behavior of the swimmer is of prime impor- 
tance in designing artificial micro-machines. Here we see 
that by changing the parameters of the system we can 
do this favor. In figure [3l we have shown that the overall 
swimming direction is sensitive to the initial phase </>o- 
Additionally and as an another example, in figured we 
have shown that by changing Lq, the length of long rod, 
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the average swimming velocity can be changed. 

As the geometry of the two-sphere swimmer is not sym- 
metric, the far field distribution of fluid velocity at the 
leading order of approximation, resembles a velocity field 
of a single force dipole. This is the main characteristic of 
most swimming microorganisms with, dipole-like velocity 
pattern. 

In summary, inspired by bacterium swimming, we 
proposed and analyzed a swimmer, constructed by two 
joint spheres. We have shown that this simple three- 
dimensional swimmer is a model for a low Reynolds pro- 
peller, captures a number of dynamical features in mi- 
croorganisms. It will be interesting to use this model 
swimmer and study many interesting problems like, the 
hydrodynamic interaction between such swimmers, the 
effects due to the confinement in the bounded fluids and 
also chemotaxis phenomena. Inspired by the colonies of 
microorganisms, we are extending our model to inves- 
tigate the hydrodynamic interaction in an ensemble of 
two-sphere swimmers. 

Appendix A: Mathematical Details 



en = (M xx - l)an + M xy a 21 + M xz a 31 + (m, - z )b 21 

~K ~ 2/0)^31 + G xx 
C12 = (M xx - l)a 12 + M xy a 22 + M xz a 32 - (m y - y )b 32 

+ G X y 

C13 = (M xx - 1)013 + M xy a 23 + M xz a 33 + (m z - z )b 23 
~\-G xz 

c 2 i = M yx au + (M yy - l)a 2 i + M yz a 31 + (m x - x )b 31 

C22 = M yx a 12 + (Myy - l)a 22 + M yz a 32 + (m x - x )b 32 
-{m z - z Q )bi 2 + Gyy 

C 23 = M yx a 13 + (Myy ~ 1 ) Cl 23 + My Z CL 33 ~ (fTl Z ~ Z )6l3 
+ Gy Z 

C31 = M zx au + M zy a 2 i + (M zz - l)a 3 i - (m x - x )b 21 
~\-G zx 

C32 = M zx a 12 + M zy a 22 + (M zz - l)a 32 + (m y - yo)bi 2 
+G zy 

C33 = M zx a 13 + M zy a 23 + (M zz - l)a 33 + (m y - yo)h 3 
-{m x - x )b 23 + G zz 



Here we present the explicit form of the matrix ele- 
ments dij, bij and Cjj which were introduced in the text: 



1 xl 

^2 



0* = - — [{! + «) + {Cr-Ct)^)] 

1 / w X 0i X j . 

B « =flji= -M (Cr " c,)( ^r ) for 1+2 

b u = 
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